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ABSTRACT 

We characterize the kinematics, morphology, stellar populations and star formation histories 
of a sample of massive compact galaxies in the nearby Universe, which might provide a closer 
look to the nature of their high redshift (z > 1.0) massive counterparts. We find that nearby 
compact massive objects show elongated morphologies and are fast rotators. New high-quality 
long-slit spectra show that they have young mean luminosity-weighted ages (< 2 Gyr) and 
solar metallicities or above ([Z/H] > 0.0). No significant stellar population gradients are found. 
The analysis of their star formation histories suggests that these objects have experienced 
recently enormous bursts which, in some cases, represent unprecedented large fractions of 
their total stellar mass. These galaxies seem to be truly unique, as they do not follow the 
characteristic kinematical and stellar population patterns of present-day massive ellipticals, 
spirals or even dwarfs. 
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1 INTRODUCTION 

During the last years, many studies have found that the most mas- 
sive (M* > 10^^ M©) galaxies were more compact in the past (a 
factor of ~4 at z > 1.5) than their equally massive lo cal counter- 
parts JDa ddi et al. 2005, Trujillo et al. 2006, Longhet ti et aLll200^ 
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Tmiillo et al. 2007u.Toft et al...2007u.Zinn et al...2007L.Cimatti ( 
2008 , Buitrag oet al.l l2008l and others). High-z massive compact 
objects are the center of many scientific debates as there is 
not a consensus in the community about so me of their proper- 
ties, l ike their ages or velocity dispersi o ns (ICenarro & Truiillol 
2009', 'van Dokku m Kriek & Franxl l2009l . ICappellari et al.l l2009l. 
Mancini et al. 20ld and others). Thus a detailed kinematical and 
stellar population analysis of these objects is crucial to under- 
stand their nature. Furthermore, it is mandatory to address the 
question of how these high-z massive compact galaxies have 
evolved into the present-day massive population. There are sev- 
eral scenarios proposed to explain this evolution, the most pop- 
ular one suggests that "minor mergers" are the dominant mecha- 
nism for the size and stellar mas s growth (Khochfar & Silkll200d, 
Naab Khochfar & Burkertll2006l. iH opkins et al. 2009, Wuvts et al. 



201QlJTmiillo. Ferreras & de La Rosall201lLlOser et aLll2012l) . This 
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scenario predicts that we should find few massive relic com- 
pact galaxies i n the local universe having old stellar populations 
(Hopkins et'aDl2009h . 

Indeed, iTmiillo et al.l (120091) (T09 hereafter), found a sample 
of 48 of these superdense massive compact ga laxies in the SPS S 
DR6 spectroscopic survey at z< 0.2 (see also ITavlor et al.ll201Qh . 
An analysis of the spectra for 29 of these objects showed that they 
had large velocity dispersions, small effective radii and young 
stellar populations. This last result is in contradiction with the 
idea that these objects are the relics from the high-z Universe. 
Unfortunately the modest quality of the SDSS spectra did not allow 
them to perform a more detailed study of their stellar populations 
(e.g, gradients, star formation histories). T09 showed the need 
to fully characterize these local compact galaxies, not only to 
understand their formation, but also to understand the formation 
and evolution of their high-z compact analogues. 

Aimed at characterizing these nearby objects, we have ob- 
tained new, high-quality spectra for seven of these compact galax- 
ies. They are representative of the T09 sample, as they cover a 
range in stellar masses, magnitudes and redshifts, so their inferred 
properties can be extended to the majority of the nearby massive 
compact objects. We performed a detailed morphological, kine- 
matical and stellar populations analysis. We derived mean SSP- 
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equivalent ages, metallicities and abundance ratios from the ab- 
sorption lines measurements and also studied their Star Formation 
Histories (SFH) using a full spectrum-fitting approach. In addi- 
tion, we obtained high sub-kpc quality imaging with adaptive op- 
tics at GEMINI-N for some of our candidates, which confirmed 
their compactness and showed that their stellar surface brightness 

t irofiles are like the ones of massive compact galaxies at z ~ 2 
Truiillo. Carra sco & Ferre-Mateu"2012'). 

The layout of this paper is the following. In Sec. 2 we introduce 
the data and the reduction process. Section 3 describes the derived 
morphologies and the stellar kinematics, while Sec. 4 focuses in the 
stellar population analysis: line strengths (Sec. 4.1), the analysis of 
the gradients (Sec. 4.2) and the star formation histories (Sec. 4.3). 
We finally discuss both the dynamical and stellar masses of these 
galaxies in Sec. 5. We have adopted the cosmology of 0^=0.3, 
Qa= 0.7 and Ho= 70 km s-^ Mpc-^ . 



2 THE DATA 

The galaxies analyzed in this article were taken from the sample of 
29 local bonafide massive compact galaxies from T09. Galaxies in 
that s ample were chosen from the NYU Valu e-Added Galaxy Cat- 
alog telanton et aPbOOllBlanton & Roweis"2007. B07 hereafter) 
with < z < 0.2, M, > 9.2 x lO^^M© and eff'ective radius < 1.5 kpc, 
adopting a Chabrier (2003) initial mass function (see T09 for more 
information about the sample selection). Note that the stellar mass, 
size and redshift were the only parameters considered for the sam- 
ple selection in T09 and that we did not apply any further cut based 
on morphology, colour, or star formation rate. We obtained high- 
quality long-slit spectra for seven massive compact galaxies in the 
4.2m William Herschel Telescope during 4 nights (22-23 July 2009 
(run A) and 3-4 June 2010 (run B)). We used the blue arm of the 
ISIS spectrograph with the grating 600B, which gives a wavelength 
coverage from 3800 to 5300 A at resolution 1.74 A for the first run 
(slit- width of 1") and 2.64 A for the second one (slit- width of 1.5"), 
with a typical seeing of 0.7''. We changed the setup on the sec- 
ond run to improve the quality of the spectra by slightly opening 
the slit. Several exposures of 30 minutes were taken (see Table 
1) for each galaxy depending mainly on their redshift, with the 
slit positioned along the major axis . Spectrophotometric standards 
and stars from the MILES library (ISanchez- Blazquez et al.l l200d 
ICenarro et al.' 2007) were observed each night for relative flux cali- 
bration. We performed a standard data reduction (bias substraction, 
flat-fielding, cosmic-ray removal, C- and S -distortion correction, 
wavelength calibr ation, sky substraction and flux calibration) us- 
ing REDUCEME (Cardiel 1999), an optimized reduction package for 
long- slit spectra that gives the error that is propagated as the data is 
handled. 



3 MORPHOLOGIES AND STELLAR KINEMATICS 

3.1 Morphologies 

We u sed the 2-D fitting code GALFIT JPeng et al.ll2002lJPeng et"aD 
12010) to obtain the eff'ective radius, the Sersic index and the axis 
ratio (see Table [T])- The program makes a convolution of the Sersic 
model with the PSF of the images, determining the best fit when 
comparing the convolved model with the galaxy surface brightness 
distribution that minimizes the x^ of the fit. The Sersic index n is 
a measure of the shape of the surface brightness profile, giving an 



idea of the morphology of the object. In the nearby Universe, galax- 
ies with n <2.5 are nor mally disk-like objects w hile n >2.5 is in- 
dicative of spheroidals ('Ravindranat h et al.1l200 2). We have scaled 
the eff'ective radius to the physical scale in kpc relative to their red- 
shift and then circularized it: 



Re - Tg ^Jb/a, 



(1) 



where R^ is the circularized eff'ective radii and r^ and b/a the ef- 
fective radii along the semi-major axis and the axial ratio derived 
from GALFIT. Figure 1 shows, for illustration, the galaxy images, 
together with the model from GALFIT and the corresponding resid- 
ual for two of the galaxies (54829 and 321479). The rest of galaxies 
can be found in Appendix A. The best fit (lowest ;^2) was visually 
inspected, to check its quality. The majority of the galaxies in our 
sample are elongated (b/a< 0.6), showing disk-like shapes. How- 
ever, two of them (galaxies 54829 and 896687) are practically cir- 
cular. In contrast, Sersic indices are mainly all above 3, suggesting 
that they are spheroid-like objects. With the resolution (~ 1 arc sec) 
and depth of the SDSS data we are not able to trace any evidence of 
merging activity, although some residuals are found in the centers. 

3.2 Stellar Kinematics 

For the kinematical study we used the Penalized Pixel Fitting 
method (PPxF; Cappellari & Emsellem 2004) with the model 
templates from Vazdekis et al.l i2010 (VIO hereafter). To measure 
the radial velocities we binned in the radial direction the spectra 
to ensure a minimum S/N(A)=10, while we binned to achieve 
S/N(A) ~ 20 for the velocity dispersion. This means that we 
summed at least two pixels for each bin, to take into account the 
eff'ect of the seeing. In those cases where the outermost radii did 
not reach the required S/N we omitted those points for the analysis. 
We find that five of our galaxies are rapidly rotating (see Fig. 
2a). The galaxies 890167, 321479 and 796740 show both the 
largest radial velocities and higher velocity dispersions (cr). For 
the remaining four galaxies we cannot find a significant rotation 
and for some of them (54829, 685469 and 2434587), we do 
not reach the flat part of the rotation curve and, therefore, we 
cannot calculate Ymax- Considering that our galaxies could have a 
significant disk component due to their low b/a, a correction for the 
inclination eff'ect is explored based on the inclination derived from 
the axial ratio. Nevertheless, we do not have enough information 
to select which galaxies are face- or- edge on. Therefore, the 
applied correction gives us an upper limit to the rotation in case a 
desk exists. As shown in Fig. 2b, the maximum velocities derived 
from the corrected curves, are larger than those derived from 
the uncorrected ones. Not reaching the Ymax for some of them 
prevented us to conclude whether these objects are also supported 
by rotation. Overall, we find high radial velocities (V^ as high 
as 200kms~^) and high velocity dispersions (~ 200km s~^) for 
most of our galaxies. The high cr is in good agreem ent with the 
values found in T09 and in lValentinuzzi et aP (l2010h for a sample 
of massive, similarly sized compact galaxies in local clusters 
(z~0.05). 

Figure 2d shows the p o sition of our objects i n the anisotropy 
diagram JBinnevI |2005|, ICappellari et al.l I2007L and references 
therein). This diagram relates the ratio between the ordered and 
random motion (Y/cr) to the galaxy observed flattening (s). The 
observed Y/cr value is obtained from the central velocity dispersion 
(equivalent to approximately 1 Rg) and Ymax- For those galaxies 
for which Ymax is not reached, we use the highest radial velocity 
obtained (Whigh), which represents a lower limit. A dichotomy 
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Table 1. Galaxy properties 



ID NYU R.A.(J200) Decl.(J200) Re(kpc) Sersic index b/a 



#exp S/N cr(kms-0 



54829 


232.57104 


-0.4885095 


1.12 


4.60 


0.90 


0.721 


0.085 


3 


41 


136.8 ±3.6 


321479 


320.21978 


11.120310 


1.20 


5.80 


0.51 


1.175 


0.128 


4 


62 


221.2 ±5.1 


685469 


335.41803 


13.987279 


1.48 


3.03 


0.45 


0.599 


0.149 


5 


70 


203.9 ±9.9 


796740 


221.90155 


43.496021 


1.24 


2.40 


0.35 


0.566 


0.182 


4 


26 


203.3 ±7.6 


890167 


234.89197 


44.297863 


0.83 


3.72 


0.63 


0.547 


0.143 


7 


32 


233.1 ±11.1 


896687 


218.94667 


54.591381 


1.63 


5.44 


0.95 


0.419 


0.130 


4 


50 


222.9 ±5.3 


2434587 


169.24737 


17.154811 


1.13 


5.45 


0.40 


0.838 


0.172 


7 


46 


205.6 ±3.1 



This table shows the main properties of the galaxies explored in this study. See T09 for other properties. The values oix^ is the reduced one from GALFIT. 



G^lfjl Model 



Residual 




Figure 1. Thumbnails for the compact galaxy 54829 (upper panels) and 321479 (lower panels). The first panel is the galaxy image from SDSS, the second one 
the fitted model from GALFIT and finally, the residuals from the fit are shown in the third panel. The white line on the upper right comer of first panels is the 
equivalent size of lOkpc at the redshift of the object. The rest of the galaxies can be found in Appendix A. 



for the elliptical family has been suggested: from one side, the 
most massive galaxies, which are slow rotators, metal-rich, with 
a flat central luminosity profile, showing evidence of triaxiality; 
and from the other side, the less massive ones, being fast rotators, 
metal-poor, with a clumpy luminosity profiles, maybe containing 
disks and being axisymmetric (e. g., IPavies et aP ll983L Ib ended 
I1988L iKormendv & Bendeil Il996h . The two lines in Fig. 2d 
correspond to models for edge-on oblate galaxies with different 
anisotropics: the upper dotted line corresponds to isotropic models 
6- 0, the lower dashed one to the linear relation 6- 0.7 Sintrinsic, 
which approximately traces the lower envelope described by 
the location of the observed fast-ro tating galaxies on the (Y/cr) 
diagram (see I Cappellari et al]|2007l for a full discussion of these 
relations). Fast rotators would lie, in this diagram, to the left of the 
dashed line. As it can be seen, the majority of the massive compact 
galaxies, without correcting them from inclination effects (big 
symbols) are fast rotators, except galaxies 685469 and 2434587 
that fall in the slow rotators regime. However, if an inclination 
correction is applied, the latter galaxies practically reach the 
separating line. Moreover, as the Ymax was not reached for some of 
them, their position should be considered as a lower limit, meaning 
that the majority of our galaxies are fast rotators. 



4 STELLAR POPULATIONS ANALYSIS 



4.1 Central Properties 



We performed a stellar population analysis using the VIO mod- 
els and the newly-defined LIS-8.4A system of indices intro- 
duced in that paper. This system is characterized for having a 
constant resolution (FWHM=8.4A) across the whole wavelength 
range and a flux-calibrated response curve. All our line- strength 
measurements, including the Lick indices ([Worthev et all 1 1994 
IWorthev & Ottavianilll99l iTrager et al.lll998k are performed on 
this system. We also used the newly-defin ed metallicity-insensitive 
index ¥1/3 o (ICervantes & VazdekisI 12009). This index, combined 
with a metallicity sensitive index, as Mgb or Fe4383, gives a very 
orthogonal model grid, given its reduced sensitivity to metallicity 
when compared to the classical HjS Lick inde x. The stellar popu - 
lation parameters were derived with rmodel ( Cardiel et al.ll2003h . 
a program that interpolates in the SSP model grids. We typically 
used the bivariated fit, although a careful visual inspection was also 
carried to check the results. 

For the central stellar population analysis we summed up the spec- 
tra within ^\Re- Before summation, the galaxy spectra were cor- 
rected pixel by pixel for the radial velocity and then all pixels were 
broadened to the central velocity dispersion. For all g alaxies, we 
corrected from emission line effects using GANDALF JSarzi et al.l 
I2OO6I) with the VIO model templates, although we only found neb- 
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Figure 2. Radial velocities derived with PPxF: first panel shows the observed radial velocities, while in the second one, radial velocities corrected from 
inclination effects are plotted. It can be seen that galaxies 896687 and 54829, which do not show strong rotation in the first panel, when corrected, they show 
slightly stronger rotations; b) Velocity dispersions from PPxF; c) The anisotropy diagram (V/cr, £). Big symbols correspond to the objects without correcting 
from inclination effects, while the small symbols are the same galaxies but correcting the inclination effect, in case they should be corrected. The dotted 
black line corresponds to the location of models for oblate edge-on isotropic {6= 0) galaxies, the dashed one to the linear relation 6= O.lstntr- Relations from 
Cappellari et al. (2007). For those galaxies in which we do not reach the flat part of the rotation curve, these values should be considered lower limits. 



ular emission in one galaxy, 685469. Figure 3 shows the index 
HjSo versus the combined metallicity indicator [MgFe50], that is in- 
sensitive to Qf-enhancement ( Kuntschner et alJl20 06). Our compact 
galaxies have young mean SSP-equivalent ages (<2Gyr). These 
results are in good agree ment with the analysis pe rformed in T09 
but in disagreement with lValentinuzzi et alJ ilOlu ), who found old 
stellar populations for their local compact massive galaxies. Over- 
all, we find mean metallicity values above solar for our compact 
galaxies, though there is a scatter in the metallicity inferred in the 
panels of Fig.4. The galaxy 685469 is not shown in these figures. 
Instead, its age and metallicity were derived from the full spectrum 
fitting with one SSP from Ulyss (see Sect 4.), due to the strong 
emission lines affecting the H/?^ and Fe5015 features. This galaxy 
is very different from the rest, being the oldest one with 2.65 Gyr 
and with a total metallicity [Z/H]= -0.55. For comparison, we have 
overplotte d some control elliptic als with similar velocity disper- 
sions from lSanchez-BlazQuez eFa l. (2006) (NGC2329, NGC3379, 
NGC4621 and NGC5812; green cros ses) and some spir a ls ob- 
served for the SAURON project in iFalcon-Barroso et al] fcOl ll) 
(NGC4235, NGC5689 and NGC6501; yellow crosses). Both con- 
trol samples are in the region of the grid corresponding to old ages. 
Thus, our objects are much younger than both ellipticals and spirals 
of the local universe. This result is also shown in Figure 5, where 
we present the inferred ages and metallicities from our index-index 
plots. Our compact galaxies (purple) show much younger SSP- 
equivalent ages and higher total metallicities than the elliptical 
(green) and spiral (yellow) control samples. 
The metallicity values derived from the different panels can be used 
to estimate, on a relative scale, the abundance ratio pattern of these 
objects. Note that the values obtaine d this way ([ZJZfe]) are a 
good proxy for the abundance ratios (lYamada et al.l i2006, VIO). 
The CN was only measured for four out of the seven galaxies due 
to the limited spectral coverage. For these four galaxies, we find 
that [ZcN/'Z^FeKO, which is at odds with the typical [CN/Fe]^0 
found in giant ellipticals of the same mass (see Fig. 4). Note, how- 
ever, that the [CN/Fe ] shows a dependence on the environment 
(ICarreteroetal.ll2004 ISanchez-Blazauez et alll2003h . with lower 
CN2 for denser environments, but our objects do not belong to any 
known cluster or group. The [ZcJ'Z.Fe] is also underabundant for 



most of our ga laxies. This result is the s ame that the one found in 
IVazdekis efaD [l991) and Cenarro et aP (^004) for massive ellipti- 
cals. This behavior is not yet clear, as the Ca, which is an of-element, 
is supposed to track Mg. For the [Zc/Zpe] and [ZMglZpe], the behav- 
ior is more complex. Visual inspection suggest that, discarding the 
three youngest objects that fall in the fully-degenerated part of the 
grid, the other galaxies do not show a significant departure from 
scaled-solar rate. However, the degeneracy of the models in these 
age regime prevents us to quantitatively determine the abundance 
ratio for these elements. We can conclude that the abundance ra- 
tio estimates for our compact objects, do not resemble those seen 
in ellipticals of similar velocity dispersion, except for the [Ca/Fe] 
abundance ratio. 



4.2 Gradients 

Taking advantage of the spatial information provided by our high- 
quality long-slit spectra we also carried out a similar analysis to 
study possible gradients in the stellar population parameters. The 
previous analysis was done using a central aperture, but we have 
extracted two other annular apertures or radii reaching, in a cou- 
ple of cases, 3 R^ and 4 R^. Table 3 specifies the radial coverage of 
each aperture, which will be used throughout the paper (rO for the 
central aperture, rl for the first annular aperture, and r2 for the out- 
ermost aperture). Overall, we do not find any strong age or metal- 
licity gradient, as shown in Figure 6. There is a trend suggesting 
the centers to be slightly younger, but these differences in age are 
really small (< 0.5 Gyr). There is also a small trend of the centers 
to have slightly higher metallicity, but these changes are within the 
errors. However, gradients in the abundance trends are difficult to 
interpret because the grids are too degenerated for such low ages, 
as stated in the previous section. 



4.3 Star Formation Histories 

Here we investigate whether the young mean SSP-equivalent ages 
obtained for these compact galaxies are due to the contribution 
of recent starbursts that mask an old embeded population, which 
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Figure 4. The age-sensitive indicator li/3o is plotted against various metallicity indices in the LIS8.4A flux-calibrated system. The SSP model grids of VIO are 
plotted. Age (in Gyr) increases from top to bottom as indicated in the labels, and metallicity from left to right ([Z/H]=-0.71, -0.40, 0.00, -1-0.22, except for the 
Ca4227 grid, where we also plot the models with [Z/H]= -1.01). Symbols correspond to Figures 2 and 3. 



Table 2. Luminosity- weighted Ages and MetaUicities 



ID NYU 



age(Gyr) [Z/H] 

[MgFe50]-I^, 



age(Gyr) [Z/H] 

Fe4383-HySo 



age(Gyr) [Z/H] 

Mgb-I^, 



age(Gyr) [Z/H] 

C4668-I^o 



age(Gyr) [Z/H] 

Full Spectral Fitting 



54829 
321479 
796740 
890167 
896687 
2434587 




1 30+0-^^ 
1 05+^-08 
1 48+^-21 

l.H-O ^ 1Q 



90+007 



35+0-05 
50+007 
20+005 



1.84+ 



50+004 

-0.07 



0.15+ 



1.30+0-06 


0.40+0-|0 


1.27+0.018 


0.17+0.004 


1-11 06 


0.35+0:10 


1.04+0.005 


0.22+0.001 


1.50+0^27 


0.10+o:|o 


1.40+0.008 


0.16+0.009 


0.80+0-08 


0.00+0:12 


0.65+0.006 


0.06+0.013 


0.99+0-05 


0.30+o:o8 


1.20+0.011 


-0.07+0.011 


1 70+032 


20+0:05 


1.60+0.009 


0.22+0.001 



Ages and metallicities derived from difl'erent line- strengths measurements (columns 2-11) for a single SSP and also inferred from the full spectral fitting 
technique of ULySS (columns 12-13). Errors were estimated with 1000 Monte-Carlo simulations using the errors on the indices and deriving Icr error 

contours in the age-metallicity space. 



Table 3. Annular Apertures coverage 



ID NYU 


rO(kpc) 


rl (kpc) 


r2(kpc) 


54829 


0.0-0.9 


0.9-1.4 


1.4-2.0 


321479 


0.0-1.1 


1.1-1.9 


1.9-2.7 


685469 


0.0-1.0 


1.0-1.8 


1.8-2.5 


796740 


0.0-1.5 


1.5-2.5 


2.5-3.4 


890167 


0.0-1.8 


1.8-3.0 


3.0-4.3 


896687 


0.0-1.1 


1.1-1.8 


1.8-2.5 


2434587 


0.0-1.2 


1.2-1.6 


1.6-2.8 



Radial coverage of each aperture. The notation will we followed through 

all the paper: rO for the central aperture, rl for the first annular radius, and 

r2 for the second and outermost aperture. 



is dominant in mass, or, alternatively, whether they are genuinely 
young objects. For this purpose we apply the full spectrum-fitting 
approach to estimate their Star Forma tion Histories (SFHs). Th e 
code used in this section is STARLIGHT JCid Fernandes et alJ2005h . 
with the VIO SSP SED library. We have defined two bins for the 
ages: "young" for lower than the SSP-equivalent age (< 2 Gyr) and 



"old" for > 5 Gyr. We have selected these two bins as the majority 
of the galaxies do not present populations in between. Table |4] (see 
also Fig.7 and Appendix B) shows the percentages derived from the 
SFHs. We can separate our objects in two families: (1) those galax- 
ies containing mainly old stellar populations (with a young contri- 
bution lower than 25% in mass); (2) those with a large fraction of 
young stars (contributing with more than 25% in mass). This last 
result is surprising, as we would not expect these massive galaxies 
to form nearly half (and even more) of their mass in recent bursts. 
The uncertainties on the L- and IVL-weighted "young" parameters 
were estimated with 2283 pairs of repeated SDSS ETG observa- 
tions. For all these duplicate spectra, we measured the parameter 
with STARLIGHT, using the same set-up and set of model templates 
as for the primary spectra. A linear relation is found between the 
Ayoung and the average S/N of each pair. The S/N of the observed 
spectra are fed into the previous relation to estimate the rms value 
of the "young" parameters presented on Table 4. 
The lack of strong gradients derived from the line- strength analy- 
sis is confirmed, as fractions remain unchanged through the diflfer- 
ent apertures. This result also implies that the starburst takes place 
along the whole structure of the galaxy. 
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Figure 3. The age-sensitive indicator HjSo. is plotted vs the composite index 
[MgFe50], which is a good proxy for the total metaUicity. All indices are 
in the LIS-8.4A flux-calibrated system and the SSP model grids of VIO are 
plotted. Age (in Gyr) increases from top to bottom as indicated in the labels, 
and metaUicity from left to right ([Z/H]= -0.71, -0.40, 0.00, -fO.22). Sym- 
bols are like in Fig. 2 and green crosses correspond to the control ellipticals 
from Sanchez-Blazquez et al. (2006, paper II), that have velocity disper- 
sions similar to those from our objects. Also overplotted as yellow crosses, 
the control spirals with simil ar velocity dispersions from the Sauron Project 
( iFalcon-Barroso etal]|201lb . Note that all the compact objects show mean 
SSP-equivalent ages smaller than 2 Gyr and they all show total metaUicity 
above solar. 
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Figure 5. Derived stellar population parameters plotted v^ the velocity dis- 
persion for our compact galaxies (purple circles). Control ellipticals (green 
diamonds) and control spirals (yellow diamonds) are also shown. Note that 
for a given velocity dispersion, our compact objects do not show similar 
ages as their local large- sized counterparts, either elliptical or spiral galax- 
ies. Moreover, their total metallicities are systematically richer. 



Apart from ST ARLIGHT, we repeat ed the analysis w ith two other 
codes, ULySS (Koleva et al' 'IOOQ") and STECKMAP (JOcvirk et all 
|2^6a, Ocvirk et al. 2006b), all with the VIO SSP SED library. The 
analysis was repeated by different coauthors independently to dou- 
ble check the possible bias introduced by the selection parameters, 
the use of different wavelength spectral regions, etc. We find that 
the three codes provide very consistent results. In order to trust the 
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Figure 6. The age-sensitive indicator Ufio is plotted against the mean metal- 
Ucity indicator [MgFe50] for our seven compact galaxies. We only show the 
model grid corresponding to the age range covered by our galaxies. Age (in 
Gyr) increases from top to bottom as indicated in the labels, and metaUicity 
from left to right ([Z/H]= -0.71, -0.40, 0.00, -fO.22). Symbols are like in the 
previous figures, but the size indicate the different apertures, with decreas- 
ing size for increasing aperture. See Table 3 for the corresponding sizes of 
each aperture on each object. 



recovered SFHs and check the reliability of the results we have 
also performed a couple of tests, described in Appendix B. The 
first test constrains the SFH derived with the full- spectrum fittings 
techniques in terms of the derived fractions. Following the results, 
we find that all our galaxies show fractions compatibles with their 
mean SSP-equivalent age (see Fig. B2). The second test checks the 
robustness of the results, as it compares the mean age and metal- 
licity from the galaxy and the one inferred when measuring the 
line- strengths of the fitted model (Fig. B3). We are able to recover 
the same age and metaUicity, although the latter with an error of 
-O.lOdex. 



5 MASS ESTIMATES 

Finally, we have revised the mass estimates for these objects, as the 
mass was one of the selection criteria of the original sample. We 
have estimated the stellar and dynamical masses following differ- 
ent approaches, as described below. The results are summarized in 
TableEl 

As mentioned before, our galaxies were selected from the NYU 
VAGC (B07) with a mean stellar mass for these seven galax- 
ies of M*507 ~9.2x10^^Mo. However, the catalog from the 
SDSS/DR7/MPA (http://www.mpa-garching.mpg.de/SDSS/DR7/ , 
hereafter MPA), gives systematically lower masses, with a mean 
stellar mass of M'^mpa ~ 8.3 X 10^^ M©. The mass estimates of 
these two analysis were obtained by fit ting the broad-band SDSS 
(u,g,r,i,z) photometry, with the models of lBruzual & Charlotl(l2003h 
(BC03, hereafter) and a Chabrier IMF (Initial Mass Function). 
However, in B07 they added the near-IR phot ometric bands of 
2MASS (J, H, Ks). Following the approach of iKauffmann et alj 
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Figure 7. Results from STARLIGHT for the central apertures of galaxies 2434587 and 54829. On the upper-left hand, the galaxy spectra (black) is fitted with 
the models from VIO (red), bottom-left panel shows the fitting residuals. The recovered Star Formation History are also shown, in light (upper-right panel) and 
in mass (down-right panel) for each galaxy. All our objects show recent bursts, some of them leading to a big contribution from the young stellar populations. 
The rest of galaxies SFHs can be found in Appendix B. 



OOOSh based on the spectral fitting, instead of the photometry, we 
used the SSP model SEDs of VIO with a Kroupa Universal IMF. 
For this purpose we use the M/L corresponding to the SED model 
that best fits the galaxy spectra, which provides a mean luminosity- 
weighted age and metallicity. The stellar mass was inferred by cal- 
culating the M/L in the SDSS r-band, adopting the GALFIT magni- 
tude as luminosity estimate, after correcting for galactic extinction. 
The stellar masses computed in this way take into account the mass 
of the luminous stars and the remnants which, for the ages we are 
considering (~ 2 Gyr), is of the order of 0.6 times the initial mass 
of gas converted into stars. The mean stellar mass calculated this 
way is M^^^p ~4.3 x 10^^ M©. In a diflferent approach, we con- 
sider instead of an equivalent-SSP stellar population, the spectra 
resulting from integrating the contribution of diflferent SSPs, with 
diflferent M/L, according to the SFH derived with STARLIGHT (see 
Section 4.3), leading to a mean W^sfh ~ 1.5 x 10^^ M©. For com- 
parison, we carried a similar approach with the SDSS fiber spec- 
tra, but with a longer wavelength-range, using the mass-weighted 
SFHs solution from STARLIGHT, obtaining a mean stellar mass of 

From the above estimates it is clear that the stellar mass de- 
termination is not straightforward. It has long been reported 
that the methodology, the models, the selection of the IMF and 
other parameters, strongly aflfect these determinations (see e.g.: 
iBell & de Jon? '200 1'. ' Kauff^mann et all [20031 . ISwindle et alJl201ll 
and Cappellari et al. 20121). For example, B07 gives a mean mass 
of M*507 ~ 9.21 X W^ Mq for a Chabrier IMF that translates to 
M*507 ~ l.Ox 10^^ Mq when employing a Kroupa, more similar 
to our M^sFH- The masses derived with the SFHs are systemat- 
ically larger than the equivalent-SSP measurements. Since in the 
SSP-method, the ages are luminosity weighted, the final values are 



strongly biased to the youngest stellar component present in the 
galaxy and, therefore, the M/L ratios are also biased towards lower 
values, as their stellar masses. The SFH derived from full-spectrum 
fitting instead, reveals an important presence in terms of mass frac- 
tions of old stellar populations, whose M/L is considerably higher 
than those of the SSP corresponding to the luminosity- weighted 
age. For this reason, the M/L based on SFH tend to be larger than 
the ones derived with a single stellar population. Indeed, the objects 
with the largest contribution from old stellar populations (321479, 
796740 and 896687 and 2434587) are the ones showing the largest 
diflferences with respect to the masses derived by means of the SSP- 
equivalent parameters (by a factor 4). Instead, the three galaxies 
with the highest fraction of young stellar population show a more 
modest diflference between the two estimates (by a factor of 2). 
We have tested the M/L obtained from the SFHs derived with two 
full- spectral-fitting codes, STARLIGTH and STECKMAP, finding con- 
sistent results. As most of the results presented here come from 
STARLIGTH, we use the masses computed from their SFHs. More- 
over, we have calculated the M/L inferred by a single SSP and by 
the SFH on our control elliptical galaxies. As their luminosity and 
mass-weighted ages are similarly old, we should not find a relevant 
diflference in the M/L, as we do. 

The dynamical mass is considered to be a good proxy for the mass 
of a galaxy. The estimate comes from the Scalar Virial Theorem 
(IBinnev & TremainduOOSL and references therein). As our objects 
show both velocity dispersion and rotation motions (the latter to a 
lesser extent), we m ay account for both contributions for the dy- 
namical mass CEpinat et aLll200^) : Mdyn^Me+Myir. These contri- 
butions are shown in Table 5. The rotationally supported term is 
the one contributing less for our galaxies: 



Me = [V^R]/G, 



(2) 
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Table 4. Starlight young population fractions from SFHs 



ID NYU 


Aperture 


%youngL 


%youngM 


54829 


rO 


71±20 


28±21 




rl 


72±19 


36±19 




r2 


73±19 


34+19 


321479 


rO 


70+18 


17+18 




rl 


68±17 


16±17 




r2 


68+17 


16±17 


685469 


rO 


82±19 


65±19 




rl 


88±19 


68±19 




r2 


89+19 


66±19 


796740 


rO 


51±19 


07±19 




rl 


52±18 


07±19 




r2 


59±15 


11±16 


890167 


rO 


99±14 


97±15 




rl 


98±10 


88±11 




r2 


97±9 


75±10 


896687 


rO 


60±18 


12±18 




rl 


55±15 


10±16 




r2 


55+14 


10+15 


2434587 


rO 


60±20 


22±20 




rl 


62±19 


23±19 




r2 


68±18 


35+19 



Fractions of young (<2 Gyr) stellar populations retrieved from the SFHs 
of the galaxies using the STARLIGHT full-spectral-fitting code. First, the 
luminosity-weighted fractions, second the mass -weighted ones. For each 
galaxy, the first line corresponds to the central aperture, the second to the 
second aperture and the third to the outermost radius considered for this 
study. 



where V is the velocity at radius R and G the universal gravita- 
tional constant. We choose V= Nmax for those galaxies for which we 
reached the maximal velocity rotation, while we take V= Nhigh for 
the rest, considering the highest radial velocity obtained. The three 
galaxies showing the highest radial velocities (from Sect. 3.2) are 
the ones showing a larger contribution from this term Following 
the approach from lBertin. Ciotti & Del Principel (l2002b . the virial 



Mw 



where 



[K{n)(r^Re]IG, 



K{n)^ 



73.32 
10.465 + {n- 0.95)^ 



+ 0.954, 



(3) 



(4) 



where n is the Sersic index, Re the circularized effective radii in 
kpc and a the velocity dispersion in km s~^ This K(n) factor takes 
into account the geometry of the mass distribution and velocity 
projection along the line of sight. Under the assumption of ho- 
mology, Cappellari et al. (2006) determined that for a de Vau- 
couleurs ri/4 profile, this factor is equal to 5. Table |5] shows the 
dynamical masses under the assumptions of either homology and 
non-homology. We obtain a mean Mdyn{nohom) = 6.89 x 10^^ M© and 
^dyn{hom) - 7.65 X 10^^ M©, both lower than the stellar masses. The 
mass estimates for the galaxy 54829 are not taken into account to 
calculate mean masses, as it is the only one having a much lower 
cr (that contributes quadratically on the equation), which would 
strongly bias the mean values to lower masses. For consistency, 
we neither use this galaxy for the stellar mean masses. 



Figure 8 shows diflferent stellar masses versus the dynamical 
mass. Filled symbols refer to the homology virial mass, open 
symbols to the non-homology one. The dashed line corresponds 
to the 1:1 relation. The first noticeable result is the difiference 
in the stellar mass depending on the methodology adopted, as 
stated above. The lowest stellar masses are obtained with the 
SSP approach (M^ssp^ yellow circles), and they increase when 
using the SFHs derived for the galaxy (M^sdss » green stars). In 
between, we have the stellar masses from the photometric fitting 
in B07 (violet diamonds). However, the most intriguing point is 
that the derived stellar masses are larger, on average, than the 
dynamical masses (except for the M*^^^). Previous studies have 
shown that dyna mical masses are on av e rage equal or larger than 
stellar masses JCappellari et all I2006L Ivan der Wei et all l2006l . 
Ivan de Sande et alj|201ll) . However, the stu dy of four compact and 
massiv e early-type galaxies at z~ 1 from iMartinez-Manso et alj 
(l201lh showed that the dynamical masses where smaller, by an 
average factor of ~6, although they could not find a plausible 
explanation. Where does the problem come from? Are the stellar 
masses wrong or are the dynamical ones? Or both? 
On one hand, regarding to the stellar mass, the discrepancies may 
be related to the large contribution of young stellar populations, 
as discussed above. The SSP-equivalent stellar masses are in 
much better agreement with the dynamical masses, as the two 
are mainly dominated by the young stellar components and the 
comparison is consistent (dynamical masses are mainly determined 
by the velocity dispersion of the galaxy, which might be biased 
by the young populations that dominate the light). Instead, the 
SFH-based stellar masses are in worse agreement and larger than 
the dynamical masses, as the old stellar components present in the 
derived SFH l ead to an incon sistent comparison. Moreover, the 
recent paper of lCappellari et al.l(i2012) has shown that variations of 
the IMF with the mass of the galaxy should be taken into account, 
being then another source in the discrepancy. On the other hand, if 
we consider that the problem is on the dynamical mass, we have 
to draw special attention to the sizes of our objects. As stated in 
IStockton. Shih & LarsonI (|2010|), it could be that we are using a 
definition for the dynamical mass that was initially calibrated for 
normally- sized objects, while our galaxies are extremely compact. 
Should we double our galaxy sizes, we would obtain higher 
dynamical masses, then we would be closer to the 1:1 relation. 
In a forthcoming paper of Trevisan et al. (in prep.), the authors 
present the M^^„-M* r elation for 40000 ETGs from the SDSS 
(iLa Barbera et al.l '2010V showing a large fraction of galaxies with 
stellar masses greater than their dynamical masses. Moreover, 
they show that those galaxies, always present small sizes, while 
they show a whole range in ages. These results point out to the 
dynamical mass to most likely be the wrong one. 



6 DISCUSSION AND CONCLUSIONS 

In this work we present the main properties of 7 of the 29 bonafide 
compact massive galaxies identified by T09 in the nearby Universe 
(z < 0.2). We find that these galaxies represent a peculiar type of 
objects, not only from a detailed analysis of their morphologies 
and stellar kinematics, but also according to their stellar popula- 
tions. The main aim of this work is to characterize these objects 
in order to understand better their formation and evolution. We 
measure large radial velocities (up to 200kms~^) and high ve- 
locity dispersions (~200 kms~^) for these objects. Morphologi- 
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Table 5. Galaxy Mass Estimates 



ID NYU 


M^ssp 


M^SFH 


M*SD55 


M*507 


M^MPA 

[x(IOIOMo)] 


M, 


Mv,„ 


Me 


Mdynih) 


Mdyn(noh) 


Model 


VIO 


VIO 


VIO 


BC03 


BC03 


VIO 


VIO 


VIO 


VIO 


VIO 


IMF 


Kroupa 


Kroupa 


Kroupa 


Chabrier 


Chabrier 


Kroupa 


Kroupa 


Kroupa 


Kroupa 


Kroupa 


type 


spectra 


spectra 


spectra 


phot 


phot 


spectra 


spectra 


spectra 


spectra 


spectra 


range(nm) 


[380-530] 


[380-530] 


[360-735] 


SDSS-F2MASS 


SDSS 


[380-530] 


[380-530] 


[380-530] 


[380-530] 


[380-530] 


54829 


4.06 


9.61 


8.22 


8.39 


7.09 


2.56 


1.91 


0.04 


2.60 


1.95 


321479 


5.01 


20.30 


16.16 


10.10 


8.19 


6.81 


3.52 


1.78 


8.59 


5.30 


685469 


4.67 


5.45 


11.10 


9.03 


7.29 


7.23 


8.54 


0.38 


7.61 


8.92 


796740 


4.45 


23.05 


13.72 


9.76 


7.65 


6.05 


8.06 


2.32 


8.37 


10.38 


890167 


2.05 


3.21 


13.92 


8.79 


12.49 


5.11 


5.11 


2.10 


7.21 


7.21 


896687 


5.11 


23.12 


11.18 


9.23 


7.39 


8.13 


5.44 


0.02 


8.15 


5.46 


2434587 


4.31 


14.42 


10.63 


8.37 


6.77 


5.52 


3.69 


0.43 


5.95 


4.12 


MEAN 


4.27 


14.89 


12.78 


9.21 


8.29 


6.47 


5.72 


1.17 


7.65 


6.89 



Stellar and dynamical mass estimates. M*ssp represents the SSP-equivalent mass estimate, while M*sfh refers to the one computed with the mix of SSPs 
from the derived SFHs of each galaxy (see sect. 4.3), for our spectra (M^sfh) and for the SDSS spectra (M*sdss )• M*507 and M*mpa are the values from 
B07 and MPA catalogs. For the dynamical masses, the virial masses are Mv^„ and My, the first accounting for the non-homology. The rotation contribution is 
stated as Mg and Mdyn(h) and Mdyninoh) are the total dynamical masses again accounting for the homology/non-homology ( Cappellari et al.l2006h . All mean 
mass values were computed without the galaxy 54829, due to its low velocity dispersion in comparison to the other galaxies. 
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Figure 8. Dynamical mass v^. stellar mass computed with diff'erent ap- 
proaches. Filled symbols state for the homology (K(n)=5), while open sym- 
bols state for the non-homology. Yellow circles are the M*55/j, violet dia- 
monds are M*507 and green stars sltqM'^sdss ■ The dashed line is the 1:1 
relation. Mainly all the stellar masses, except for the M*ssp, are similar or 
bigger than the dynamical mass. 



cally, they are very compact (Re ~ 1.3 kpc) and with Sersic in- 
dices typically higher than 3, indicative of spheroids. However, 
according to their axial ratio, they look elongated on the sky 
(b/a< 0.6). From an analysis of their stellar populations, we derive 
young mean SSP-equivalent ages (< 2 Gyr) and metallicities solar 
or higher ([Z/H] ^ 0). The full-spectrum-fitting derived SFHs reveal 



that these objects contain an old embeded population but that they 
all experienced recent bursts of star formation during the last 2 Gyr. 
The latter contribute strongly to the total light (> 50%) and, in some 
galaxies, even to the total mass (> 30%). We divide our sample in 
two classes depending on the strength of the recent bursts: (1) Four 
galaxies show mainly old stellar populations (^75% of the total 
stellar mass); (2) and three objects contain a very important contri- 
bution of young stellar populations (> 30%). 

Although local compact galaxies show velocity dispersions 
and Sersic indices similar to those in local normally- sized massive 
ellipticals, they show important rotation curves and seem to be 
rotationally supported. IMoreover, unlike ellipticals of similar 
mass, our compact galaxy sample shows young SSP-equivalent 
ages (< 2 Gyr) and higher metallicities. 

Despite the fact that our galaxies display disky morphologies and 
strong rotation curves, they dififer from spiral galaxies in their 
Sersic indices and their velocity dispersions (the latter showing 
lower values for both; n < 2.5 and cr < 150 km s"^. IMoreover, our 
objects are younger than spirals of similar mass. 
For completeness, we compare our objects with dwarf ellipticals, 
although their masses are much lower than the ones measured in 
our sample. IVLany dwarf galaxies show young mean luminosity- 
weighted ages, which can be similar to those inferred for our 
galaxy sample, but in general they have lower metallicities. 
Kinematically, our galaxies are rotating as dwarf ellipticals, but the 
latter have lower Sersic indices and lower velocity dispersions. 
Our objects present op t ical spectra resembling 'E+A galaxies 
toressler & GunnI ll983L ICouch & SharplesI Il987h . E-hA objects 
are considered post-starburst galaxies due to the presence of strong 
B aimer absorption lines, indicative of a young population (<1 Gyr) 
and an absence of major emission lines ([Oil] or Ha) indicating 
the cessation of star formation. It is commonly accepted that E-hA 
represent the transition phase from a disky, rotat ionally supported 
system to a spheroid , pressure siipporte d one (JYang et al.l |2004 
IZabludofif et al.l 1 19961 . IPracv et al.l l2009h , showing a widespread 
of properties (ages, metallicities, velocity dispersions, sizes). 
Because of their optical spectra and their SFHs, our objects are 
post-startburst galaxies. However, it is worth noticing that none 
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of our galaxies, which we selected from the SDSS d atabase, was 
previously identified as E+ A (e.g., see the catalogs of Ib lake et alj 
l2004i lGotdl2005l . lGotd l2007^. The diff^erent catalogs take diff^erent 
limiting values on the equivalent width of the lines that determine 
the family, reason why our objects may not be found in them. 

The goal of this work was to characterize a sample of local 
compact massive galaxies, which were initially thought to be the 
relics of massive compact objects at high-z. Instead, we have 
shown here that they seem to be exact copies of those high-z galax- 
ies, as their inferred properties (masses, sizes, luminosity- weighted 
ages) are nearly the same. We must note, however, that their SFHs 
do not agree with the massive compact galaxies at high-z, which 
are genuinely young at those redshifts. Although our objects 
have experienced important recent starbursts, some of them are 
compatible with having very old stellar populations. We compared 
the properties of our sample with those of ellipticals and spirals in 
the local Universe. While they share a number of similarities, the 
compact galaxies also present many important diff'erences (e.g., 
their young ages) and they deviate from various relations expected 
in the local Universe. These diff'erences do not allow us to classify 
them neither as ellipticals nor spirals. Instead, their optical spectra 
resembles those of the E+A galaxies. May it be, that we are in 
fact seeing an intermediate stage on the thought transformation 
from spirals to ellipticals via E+As? We have found a new class of 
objects, unique and rare in the nearby Universe. 
Based on the derived SFHs, we can divide our objects in two 
groups. The first one is composed by those compact massive galax- 
ies with large fractions of old stellar populations. These objects 
could be considered the relics of the high-z universe, although 
their SSP-ages are still young, due to the non-negligible fraction 
of young components. In the second group, we include those 
galaxies with a huge fraction in mass of young stellar populations, 
analogous to the high-z massive galaxies. This group is the most 
intriguing. How is it possible to form half of the mass of a massive 
galaxy so recently (< 2 Gyr)? Whatever mechanism responsible for 
this, demands huge quantities of gas to create this amount of stars 
in such a short time. Is this mechanism analogous to that forming 
the massive galaxies at high-z (cooling flows or the merger of two 
"gas rich" spirals, e.g. Ricciardelli et al. 2010)? The absence of 
clear stellar gradients along the nearby compact galaxies suggests 
that the mechanism triggering the starburst in these objects is 
likely a global one. This could indicate that the structure of the 
nearby compact massive galaxies has been assembled during this 
major event, being the old stellar populations the stars associated 
to the merging units and the young component a product of the gas 
collapse into stars. 

Although we have characterized their stellar populations and kine- 
matics, we still need to understand their formation and evolution 
and assess the issue of why these local massive galaxies were born 
with such a compact structure. A first approach to solve these open 
questions would be to expand our wavelength range as, for the 
moment, only optical data is available. Moreover, we will have to 
wait until high-quality spectra for the high-z analogs are obtained, 
to compare their relevant stellar properties (e.g., the abundance 
ratios and total metallicity). We also require high quality imaging 
to properly loo k to their m orphologies and trace their merging 
activity (see e.g. lTruiillo. Ca rrasco & Ferre-Mateu 2012). 
The local compact galaxies analyzed in this paper seem to be 
unique and scarce in the nearby Universe and, because of their 
relative vicinity, they are a perfect laboratory to explore with un- 
precedented detail the formation mechanisms of massive galaxies 



at high-z, particularly those involving huge amount of gas leading 
to the enormous formation of stars we have measured. 
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APPENDIX A: GALFIT 

In this section we show the rest of our galaxies not shown in Section 
3. First panel corresponds to the original image from the SDSS 
database, the second is the model created with the PSF of a suitable 
star close to our object, and the third panel shows the residuals from 
the substraction. We chose these as the best results by their reduced 
X^ and by visual inspection. It is clear that for some galaxies there is 
a remaining residual in the center, which could not be removed with 
a simple Sersic profile. Bulge-disk decompositions could be made, 
but SDSS quality is not optimal, high-quality images are needed to 
properly study the morphologies of these objects. 



APPENDIX B: SFHS 

In this Appendix we show the SFHs recovered with the code 
STARLIGHT. Each full- spectrum-fitting code gives us diff'erent in- 
formation, so we choose the most suitable in each case. For exam- 
ple, ULySS is very efficient giving 1-SSP fits, as shown in Table 2. 
However, we find more useful the output from STARLIGHT when 
studying the contribution of the main stellar populations. Finally, 
STECKMAP was used as double-check for our results. All three codes 
were fed with the VIO models. For each galaxy, we have analyzed 
the SFH for the diff'erent apertures, to study possible gradients. In 
Fig. Bl, each row on the panel corresponds to a single galaxy, with 
the first column being the center (rO), the middle one the first aper- 
ture (rl), and the third one the outermost aperture (r2). Each panel 
contains the galaxy spectra (in black), the fitted spectra from the 
model templates (red) and the residuals, showing, if any, the re- 
jected pixels and masked regions (green). The SFHs are shown on 
the right of each spectra, both in luminosity (up) and mass (down). 

We have also tested the reliability of the full- spectral-fitting 
techniques to estimate ages and metallicities. We have applied two 
diff'erent tests, explained below. 

TEST 1: This test constrains the fractions of young and 
old populations recovered with the full- spectral-fitting technique. 
Given the mean SSP-equivalent age and metallicity for a galaxy. 
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Figure Al. Thumbnails for each compact massive galaxy showing for each object, the galaxy image from SDSS (first panel), the fitted model with GALFIT 
(second panel) and the residuals from the fit (third panel). The white line on the upper right comer of all the first panels is the equivalent size of lOkpc at the 
redshift of the object. 
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Figure Bl. Results from STARLIGHT for the three apertures in each galaxy (rO, rl and r2). On the upper-left hand, the galaxy spectra (black) is fitted with the 
models from VIO (red). The bottom-left panel shows the fitting residuals and the masked or rejected pixels (green). The right column on each panel shows the 
Star Formation History recovered in light (up) and in mass (down). It can be clearly seen that all our objects show recent starbursts, some of them leading to a 
big contribution of young populations. 



derived from the SSPs (Single Stellar Population) model grids, we 
can also recover the same mean age with a combination of two 
SSPs, one younger and the other older than the SSP-age. These 
two SSPs can be weighted by different proportions, either in mass 
or in light. We created a sample of synthetic spectra with exactly the 
same configuration as for our galaxies, using the MILES webpage 
tool ''Get spectra from a given SFH'\ where our input were dif- 
ferent SFHs, one containing a young population (< 2 Gyr) and the 
other containing an old population (> 5 Gyr), with varying propor- 
tions. We slightly vary the imposed limit (2, 1.5, 1 and 0.8 Gyr) to 
match the estimated mean luminosity- weighted age for each object 



in our sample. We then measure their line- strengths in the same way 
we did for our galaxies and plotted them in the Hfio vs [MgFe50] 
diagram (Figure B2a). The diflferent combinations of ages with a 
fixed metallicity are plotted in dififerent colours for [Z/H]= -1-0.22 
and [Z/H]=0.0. Only those points (in black) that lie close to our 
galaxies are the possible SSP combinations that will give the galaxy 
mean age. However, the dependence with the metallicity has also 
to be taken into account because of the known age-metallicity de- 
generacy (a young metal-rich system has a similar integrated en- 
ergy distribution (SED) as an older metal-poor one). To illustrate 
this, we plot the extreme case where a young metal-rich population 
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([Z/H]=+0.22) is combined with an old metal-poor one ([Z/H]=- 
2.43). All the gray points correspond to those possible combina- 
tions of SSPs with different ages and different metallicities that lead 
to our galaxy values. As we know the amount of old population for 
each of them, we are able to impose limits on the fractions required 
to achieve our galaxies values. At a fixed age, less metal-rich com- 
ponents contain smaller fraction of old populations. At fixed metal- 
licity, as the mean luminosity-weighted age increases, the fraction 
of old population can reach any value. For example, for the galaxy 
2434587 with mean SSP-age 2.24 Gyr, the amount of old popu- 
lation varies from 0.05% to 95%, while for the galaxies with less 
than 1.0 Gyr, only SSPs with more than 60% of old population are 
feasible. The fractions listed in Table 4 are in agreement with this 
test, except for galaxy the 54829, whose young population is higher 
than expected. 

TEST 2: This test is a simple reanalysis of the output of 
the full-spectral-fitting. We measured, in exactly the same way 
that we did for our galaxies, the absorption line- strengths of the 
best model fitted to each galaxy with STARLIGHT. We plotted 
these values in the index-index diagrams to compare them with 
the observed value from the galaxy spectra (for both our galaxies 
and for the elliptical control sample). It appears (Fig. B2a) that 
ages are well recovered and metallicities tend to be slightly lower 
by -O.IOZ©. For old stellar populations, the method based on 
the B aimer indicator vs. the metallicity indicator t end to provide 
older ages than the ful l- spectral-fitting approach (IVazdekis et al.l 
1200 iL iMendeL Procto7& F orbes 2007). This effect is seen for 
the control elliptical galaxies on the bottom panel of the figure. 
Instead, in stellar populations with with strong contribution of 
young components, the SSP-equivalent age is strongly biased 
towar ds the age of these young components (e.g.. ISerra & Tragej 
l2007h . 
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Figure B2. a. Left panel: The age-sensitive indicator H/?o is plotted against the composite metallicity index [MgFe50], to show the results of combining two 
different SSPs (young+old) with a fixed metalHcity of [Z/H]= 0.0 and 0.22 (orange/brown symbols). We also show the eff'ect of combining an old metal-poor 
with a young metal-rich SSP (horizontal sequence of faint points). In gray, we highUght the possible combinations that lead to the observed values, b. Right 
panel: I^o-[MgFe50] plots for the second test. Each panel refers to each galaxy of our sample. Large filled symbols are our compact galaxies as in Figure 
2, while the crosses refer to the values derived from the fitted spectra from STARLIGHT. For comparison we show the results for the control ellipticals in the 
bottom panel. 



